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1 Introduction 



A calculating of amplitudes of QCD high-energy scattering processes is a complex task based on 
different methods and approaches of high-energy physics, see [IJ El El IH E] • In general this task is very 
complicated, and, therefore, the simplified models which have dynamics similar to QCD ones could 
serve as a good polygon for the different ideas and methods check. This is the reason why RFT-0 
(Reggeon Field Theory in zero transverse dimensions) model is attracted much interest during few 
last years, see [JJ. The approach, which was formulated and studied a long time ago, even before 
the QCD era ( see [8]), has a dynamics which is very similar to the solutions obtained in framework 
of interacting QCD Pomerons, see [91 110] . 

In the paper |10j was shown, that we could understand the real RFT dynamics solving equations 
of motion of RFT-0. Quantum amplitudes in RFT-0 model framework are also calculable, see [7], 
and this is the only source for understanding, if only partially, of quantum dynamics of RFT based 
on QCD Pomerons. Nevertheless, even at classical level, there are some open questions which are 
interesting from point of view of QCD high-energy description of the process. These questions are 
about the asymptotic solutions of equations of motion, their stability and properties, see more details 
in [TO]. In contrast to what we have in RFT in two transverse dimensions, the classical solutions of 
RFT-0 could be investigated analytically and achieved results could be applicable also in high-energy 



The paper is organized as follows. In the next section we consider a way to approach asymptotic 
solutions of the model whereas in the Section 3 we analyze steady states of equations of motion. In 
Section 4 we discuss the interconnections between the steady-states and asymptotic solutions, relating 
the stable points of the equations with corresponding asymptotic solutions. Section 5 is dedicated to 
additional analysis of a stability of some asymptotic solutions. In Section 6 we investigate a possibility 
of periodic limited cycles of the solutions and Section 7 is a discussion of results and conclusion of the 
paper. 

2 Asymptotic behavior of equation of motion 

The Hamiltonian of our problem and corresponding equation of motions have the following form: 



QCD. 



H = fiqp — 7p 2 q — -ypq 2 



(1) 



and 



<1 = 



Hq - 7<? 2 - 2^pq 



p = -\ip + 7p 2 + 2jpq 



(3) 
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see [8], where /x is an intercept of bare Pomeron and A is a vertex of triple Pomeron interactions. 
The "time" variable of the equations is a rapidity , i.e. q = ^ with y = ln(s/so) and with s as a 
squared total energy of the process.. The boundary conditions are symmetric 

q(y = 0)=A, p(y = Y)=A (4) 

where Y is the end of the "time" interval, i.e. Y is a total rapidity of the scattering process. We 
also request, that an amplitude of the "scattering" process, which depends on q and p variables, will 
preserve the so called "target-projectile" symmetry 

Ampl.{q(y, A),p{y, A)) = Ampl.(q{Y - y, A),p(Y - y, A)) (5) 

see in [10] the definition of the scattering amplitude in the RFT. 

In order to qualitatively understand behavior of the equations at large rapidities we rescale the 
rapidity and vertex of the problem: 

y ^ fiy, X = j/ fi (6) 

obtaining new equation of motion 

q = q-Xq 2 -2Xpq (7) 
p=—p + Xp 2 + 2Xpq. (8) 
We see, that in the rescaled equation the asymptotic limit for the rescaled total rapidity Y 



could be achieved also by the limit 



Y ReS caled = V Y -> OO 



\i — > oo , 7 a ^ oo 



where also 

X = ~/ / fi oc a s — > oo . 

It means, that asymptotic solution of equations of motion for the case Y — > oo are the same as 
solutions in the limit A — > oo . Moreover, because our function must be analytical function of a s , the 
solutions will be valid also when Y — > oo at finite A. We note also, that due the initial condition 
Eq. (jU), the asymptotic behavior of the q = q(y,A) is achieved in the limit y — > Y — > oo whereas 
asymptotic of p = p(y,A) function is lies in the limit y — > with Y — > oo. 

The rescaled Hamiltonian of the problem now has the following form: 

H = qp - Xp 2 q - Xpq 2 = E (9) 

with equation of motion 

q = q - \ q 2 - 2Xpq (10) 
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p = -p + Xp 2 + 2Xpq (11) 

where 

y -»• ny E E / fi 

in equations of motion and Hamiltonian 0. The solution of the equation Eq. Q for variables q and p 
at the limit of large A are two pairs of the functions: 

( 1 E E \ 

Ql > 2 = [\- p - AV' W) (12) 

and 

/ 1 E E \ 

^ = U" g "X?' A?J" (13) 

Further, as any pair of asymptotic solution {q,p) i of equation motion Eq. (fTU|) - Eq. (fTTj) . we will 

consider the pairs Pi) i = 1..2 which contain functions Eq. (fT2|) - Eq. (fT3|) , 



3 Steady states of equations of motion 

The definition of steady states of the equation of motion is a standard: 

q (1 - \ q - 2\p) = (14) 

- p (1 - Xp + 2\q) = (15) 

Due the fact of absence of rapidity dependence in these equations we assume that at least part 
of solutions (steady states) of these equations describe the solutions of equation of motion in the 
asymptotic limit Y — > oo, whereas y — > Y for the q = q(y, A) and y — > for the p = p(y, A) 

The non-trivial steady states of the equations could be easily found 

(QSt^St^ = (p0) (16) 

(qst,PSt) 2 = (o, (17) 

( qS t,PSt) 3 = (^.^) (18) 

In general, the linearized analysis of the equations of motion around these states could be performed. 
The Jacobian matrix 

/ 1 - 2Xq - 2Xp -2Xq 
J = 19 
2Xp -1 + 2Xq + 2Xp 



1 Further, we shell use everywhere the sign E instead E / fj,. 
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lead to the following eigenvalues for these three steady states 

A{ 2 = ± 1 \f = ± 1 A* 2 



= ± — 



(20) 



Unfortunately, the fact that, for example, the first two steady states are saddles does not say a 
lot about particular corresponding asymptotic solutions which also must satisfy boundary conditions 
Eq. Therefore, in the next two section we will study these steady states in correspondence with the 
solutions Eq. (|12p - Eq. (|13p in the sense of convergence of asymptotic solutions to the corresponded 
steady state in the limit of asymptotically large rapidity. In the case of the convergence existing we 
could call a steady state clS ct stable point. 



4 The (qst,Pst)i an d {qst,Pst)2 steady states and corresponding so- 
lutions of equation of motion 



In this section we compare the solutions Eq. (|12p - Eq. (|13p with the results of the steady states 
analysis Eq. ([To]) - Eq. (fTT]) . We note, that the energy, which corresponds to these steady states is 
zero: 

E = (pq-Xpq 2 ) =0. (21) 

v ' {{qst,pst) 1 ,{qst,p.st) 2 } 

Therefore, in the limit E — > a correspondence between these states and the asymptotic solution is 
established easily: 

(qst,PSt)x = Q, 0^ = (qi,p 2 ) E =o ( 22 ) 

and correspondingly 

(qst,PSt) 2 = = (02,Pi) E =o - (2 3 ) 

In order to find the asymptotic solution which corresponds to the steady states/stable points of 
Eq. (|22p (we will look for the solution for the states of Eq. (|22p . the solution for Eq. (|23p could be 
found similarly) we will consider the equation Eq. ([7]) as the equation consisting of two parts. The 
first part provides a fully integrable equation, whereas the second part, 2 A qp , is a small perturbation 
around p oc 0. Keeping only leading terms in both Eq. ([7]) - Eq. ([5D we arrive, therefore, to the 
following system of equations: 

q = q - Xq 2 (24) 

p = -p + 2Xpq. (25) 

The consistency of the approximation we check thereafter by taking asymptotic limits in the pair (q, p) 
of solution ans comparison of obtained result with conditions Eq. (|22p or Eq. f|23|) . 
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Integration of the system Eq. (|24p - Eq. (|25p with given initial conditions Eq. @ leads to the 
following functions: 

"'^ = iA("l)4 1 (26) 

In order to find the asymptotic behavior of the solutions we take the limit Y — >■ oo and finite A. We 
see, that in this limit, where y — > Y — > oo for q and y — > , Y — > oo for p, we obtain 

These asymptotic values coincide with the asymptotic values ( <7i , P2 ) from Eq. (fT2|) - Eq. (fT3|) and 
stable points Eq. (|16p . As expected, taking firstly the limit A — > oo, we will come to the same stable 
point Eq. (g2D - Eq. <II3|I : 

, . ^4e y A->oo ^ y-^oo 1 

g(y ' A) = AA(e»-l) + l ~> A (30) 

P(V,A) = Aj-y ( j^'-lVl ) 2 X ^^°° Ae~ Y - (31) 
in the limit Y Rescaled oc a s Y -> oo. 

The second solution, which corresponds to the stable point {qst ,PSt)2 could be found similarly. 
It has the following form 

f AX ( e Y -y - l) + l\ 2 

q(y,A) = Ae y \ - V v / (32) 

yvy 1 I (e y - 1) + 1 I v ' 

p{y > A) = A\( e y-v -i) + i- (33) 

We also see, that an amplitude consisting of the pairs of solutions, Eq. (|2U|) - Eq. (f27|) and Eq. 
Eq. (|33p . will preserve the symmetry of the problem Eq. ([5]) 



5 The ( qst , Pst )s steady state and corresponding solution of equation 
of motion 

We search an asymptotic solution which characterized by following steady state: 

( qS t,Pst) 3 = (^.^) • (34) 
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As a pair of asymptotic solutions (q , p), which could corresponds to this state, we take the following 
combination 

( 1 E 1 E \ 

(ft,Pl) = [j -P- ^,~ x -q- j^) , (35) 

and our further task is a check of the self-consistency of this assumption. The Hamiltonian for this 
steady state/stable point is not zero: 

H(q=l/3X,p = l/3X) = E = (36) 

so we need to find solutions of equations of motion which provide this non-zero energy. 

Let's consider full equation for the q variable and, as we did previously above, we will expand the 
last " perturbative" term in the equation taking there p = pi from Eq. 



q = q — \q 2 — 2Xpq = q — Xq 2 — 2Xp\q (37) 

We obtain 

q = q - Xq 2 - 2X Pl q = q - X q 2 - 2 A q Q - q - -^j =-q + Xq 2 + (38) 

Again, considering the last term as an perturbation, we take there q = qst from Eq. (|34p and rescaling 
rapidity fiy — > \i X y we get finally: 

q = q 2 -q/X + 6E. (39) 
The solution of this equation with the initial condition given by Eq. Q is 

1 , A - 1/(2A) + E + e 2 ^y (A - 1/(2A) - e) 
„ _|_ \ /_ (40) 

~ 2A A — 1/(2A) + E - e^y (A - 1/(2A) - E^ 
where 



E = \llx^- 6E ^ 

and rapidity is 7 y . In the asymptotic limit y — > Y — > 00 we obtain a simple expression 



« = k- E - (42) 



For the p variable we similarly obtain 



1 , A - 1/(2A) + E + e 2ElY -y^ ( A - 1/(2A) - E 

P=^ + E ; ; ) (43) 

2A A - 1/(2A) + E - e 2E ( Y ~y) (A - 1/(2A) - Ej 
that in the asymptotic limit of variable p, which is y — > , Y — > 00, gives 

p = Yx~ e - ^ 



How, for the verification of the self-consistency of the solution, we back to the asymptotic pair of 
solutions Eq. (f35j) . We have 

1 . 1 E 
q= 2X- E =X- p -X^- (45) 
Taking value of p from Eq. (|44p and value of E from Eq. (|4T|) we obtain 

-2XE (l/(2A) - eY = E 2 /6 - 1/(24A 2 ). (46) 



Among three solutions of this equation there is a solution 

1 

6A 



E = j- (47) 



which gives 

= (s-^s (n-n) =(»■»■). < 48 > 

and 

£ = 1/(24 A 2 ) - £ 2 /6 = = if( g = 1/3A , p = 1/3A ) . (49) 

Clearly, the required amplitude's symmetry Eq. © for the Eq. (|4"0j). Eq. (|4"3"]l functions is also 
preserved. 

Thereby we proved that our solution is self-consistent ant that solution Eq. (|35p (i.e. functions 
Eq. (|4(jp . Eq. (|43p ) is related with the Eq. (|18|) stable point. Other pairs of asymptotic solutions for 
E 7^ 0, such as ( qi , p2 ) , ( <7i , P2 ) and ( (72 , Pi ) , do not satisfy the condition of correspondence of 
their large rapidity limits with the steady state Eq. ([TBI. 



6 Stability analysis and Lyapunov function 

Obtained asymptotic solutions are stable in the sense of their behavior at asymptotically large rapidity. 
Nevertheless, it is interesting to check the asymptotic stability of the solution from the point view of 
Lyapunov function, namely, we will try to construct the local Lyapunov functions for each stable point 
and corresponding asymptotic solution obtaining the stability criteria from this side of the problem. 

We will look for the local Lyapunov functions, it means that we will search for the functions which 
have following properties: 

1. first property: 

F(q,p) > (50) 
for the q and p defined in neighborhood region of stable point; 

2. second property: 

F(q,p) q=qsp , p=psp = 0; (51) 
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3. third property: 

( dF(q,p) 
V dy 

in neighborhood region of stable point. 



< (52) 



6.1 Lyapunov function for Eq. (1161) - Eq. (1171) stable points 

We consider stable points Eq. (fTU]) 

(qst,Pst) 1 = (~, o) (53) 

and corresponding asymptotic solution (qi,P2) defined by Eq. (j2~6|) - Eq. ([27]) . For the stable point 
( qst ,PSt )i we consider a following function as a candidate for the Lyapunov function: 

F(q,p) = H - p n + Xp 2 q (54) 

where H is the Hamiltonian Eq. ([9]) and n > 3. Now we will check properties of this function as a 
Lypunov one. 

1. We look for the behavior of the Lyapunov function around the stable point and we have at 
leading order at asymptotically large Y : 

F(q,p) iqi , P2 ) = (qp ~ Xpq 2 ~ p n ) ^ = ^ ~ XA) > 0, (55) 

the condition n > 3 is used here for the derivation of the answer and where we assume A A < 1 . 

2. Checking second request for the Lyapunov function we obtain 

F(q=j,p = 0) = (56) 
simply by definition of the function which is proportional to p. 

3. Third property of our function is the following 

f^M) =(qp- X P q 2 -p n )' =(H -np n - 1 p) ( (57) 

Here Hq = qp — Xpq 2 is a "fan" Hamiltonian with precise solution of the corresponding 
equations of motion in the form functions Eq. (|26p - Eq. (|27p . Therefore, we have Hq = and 
obtain: 

= _( np n~l-\ np n { _ 1 + Xp+2Xq) (5g) 

In leading order, around (q\ = j,P2 = e), we obtain: 

f^l) = _„,. = _„^ <0 . (59) 

where the value of e is taken from Eq. (|29|) . Thereby we see, that function Eq. (|54|) is the Lyapunov 
function for the stable point Eq. ()53|) . A Lyapunov function for the stable point Eq. (|17|) could be 
constructed similarly. 
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6.2 Lyapunov function for Eq. ( 1181) stable point 

For the stable points Eq. (fT8|) 



( qS t,Psth= (j-yj~ x ) (60) 

and corresponding asymptotic solution Eq. (|4U[) and Eq. (|43|) as a Lypunov we consider the following 
function 

F(q,p) = In (H/(Xq 2 p)) (61) 
where H is the Hamiltonian of Eq. ([9]). This function has following properties. 

1. In the neighborhood of the stable point, where q = ^ — e and p = <^ — e with e < 1, see 
Eq. (|4U|) and Eq. ([43 1) . we have in leading order on e: 



F(g,p) = ln(l + 9eA) ps 9eA > (62) 

with 

e 

~3A A — 2/3A 



e~ y /3 A - 1/3A 
e = ^77"^ oToT ■ ( 63 ) 



It is interesting to note also, that for 

1/3 < ^A < 2/3 (64) 
condition Eq. ([62]) for the function Eq. (I5T| is not satisfied. 



2. Taking precise value of the stable point Eq. (fTHj) we see that 



F(g=i )F =i) = lnl = 0. (65) 



3. Taking the derivative of this function over rapidity we obtain: 



V dy J q p 

At leading order on e we have that H = and we obtain 
fdF(q,p) 



(dF{q,p) \ _ k _ 2 i_P_ (66) 



V dy 

At leading order we, therefore, have 
fdF(q,p) 



2 (1 - Xq - 2\p) - (- 1 + Xp + 2Xq) . (67) 



V dy J =(-l + 3A 9 ), = £_ e =-3 £ <0. (68) 



Thereby we see, that our function Eq. (|61j) is a Lyapunov function for the stable point Eq. (|18|) . 
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7 Periodic limited cycles of solution 



In order to show the presence of periodic cycles in asymptotic solutions of our equations, we back to 
the original Hamiltonian Eq. (pQ) 

H = fiqp — jp 2 q — jpq 2 (69) 
and determines the vertices of equation correspondingly to the QCD Pomeron vertices : 

H = a s , 7 = a 2 s , (70) 
where a s is QCD coupling constant, see [6]. Introducing new rescaled fields 

p -> ^ , ± , (71) 

a s a s 

we obtain rescaled Hamiltonian: 

H = — (qp - p 2 q - pq 2 ) . (72) 

a s v ; 

Now we follow the results of article [11] . Consider the system of three variables q,p and z: 



Q = 


q (1- 


q - 


ap 


- (3z 


p = 


p(l- 


p - 


Pq 


— a z 


z = 


z(l- 


z — 


a q 


- Pv 



(73) 

with additional condition a + (3 = 2. The dynamics of following function 

F = p + q + z (74) 

could be described with the use of equations Eq. (f73|) : 

f g=F(1 - F) . (75) 

At Y — > oo the solution of this equation is simply 

F -> 1 , p + q + z ->• 1 . (76) 

Therefore, at asymptotically large rapidity Y possible solutions of Eq. (|73p lie on the triangle Eq. (|76p . 
i.e. all orbits at large enough rapidity set up on this triangle. Further we project this motion on the 
q — p plane by use of projection 

z = 1 — q — p (77) 



which has the place at large rapidity in system Eq. (|73p . Rewriting the equations Eq. (|73p we obtain 
for q — p plane dynamics: 

a ~ , 

q = ^r^q(l - q -2p) (78) 
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a - 

P = 2"^p(l-p-2g). (79) 

We see, that these equations could be considered as equations of motion for the following Hamiltonian: 

R = a_l^ qp _ p2q _ pq2 y (8Q) 

Comparing Eq. (|8Up and Eq. (|69|) we conclude, therefore, that our Hamiltonian of interests could be 
considered as projection of dynamics of the system of three variables Eq. ()73p on the plane q — p in 
the limit of large rapidity Y with the 

a = 1 + — (81) 
u s 

= 1 - — B (82) 

Us 

We conclude, therefore, that at large rapidity limit the dynamics of the system of two variables Eq. ((7J) 
and system of three variables Eq. (|73p is the same. The systems are dual in the limit of large Y. 

In order to illustrate the periodic solutions in Eq. I[73|) . and therefore periodic solutions in Eq. ([7J), 
we introduce an additional function: 

W = qpz. (83) 
The dynamics of this function is described by following equation: 

4- HW) = 3 - 3W = 3 — H F ) , (84) 
ay dy 



w{y)=wJ^X v ^wJ-±-) =C . (85) 



that gives after the integration: 

3 

\F(0) ) ' " u \F(Q) 
We see, that overall dynamics at large Y is restricted by the following conditions 

p + q + z = 1 (86) 

and 

pqz = C . (87) 

Conditions imposed on the solutions of the Eq. (|73|) at large rapidity result as cycle orbits, which are 
intersection of the triangle Eq. ([86]) and hyperboloid Eq. ([87]) . In this case, the projection Eq. ([77]) . 
determines also cycle dynamics on the plane defined by equations of motion Eq. (|7|). 

The duality of the systems Eq. (UJ and Eq. ([73]) could be also understood from the point of 
conditions needed to apply on the systems in order to obtain the solution. In the case of system 
Eq. (UJ we need the boundary conditions Eq. and value of the energy of the system E. In the 



2 It is interesting to note that the Eq. (|73f) system is drastically simplified in the limit 
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case of equation Eq. ([73]) we need only boundary (initial) conditions (qo , po , Zq , ). Excluding z from 
Eq. dTTJ) and Eq. ([57]) we obtain: 



pq (1 - q - p) = Cq = E . 
We see that the following relation for the dual systems exists: 

E=C = q0P0Z ° . (89) 

qo + po + z 

Hereby we conclude, that the energy parameter in system Eq. (JTJ is determined by the initial value 
zq in the system Eq. ([73]) at given (go , Po)- 

In our particular example for stable points Eq. ([22]) at E = we see, that po(y = , Y — > oo) — > 
provides E = Co = for corresponding asymptotic solution Eq. ([26]) - Eq. ([27]) through Eq. ([55]) . 
In this case, in order to satisfy Eq. ([86]) - Eq. ([57]) we also need to conclude that z{Y — > oo) — > . In 
general, therefore, whole dynamics of the system Eq. ([7J at large rapidity and at E = is described 
by the line 

q(y) + p(v) = i (90) 

with stable points given by Eq. (j!6|) - Eq. (|17|) . In the case when E ^ at the limit of large Y, the 
cycles are not lines anymore but they are cycles determined by Eq. ([86]) - Eq. (jSj 



8 Conclusion 



In this paper we investigated the form and behavior of solutions of RFT-0 theory at the limit of 
asymptotically large rapidities. We developed a method which allows to find an analytical form of 
asymptotic solution basing on the analysis of steady states of the equations of motion of corresponding 
Hamiltonian. Found solutions satisfy the boundary conditions of the problem, Eq. ([!]), and preserve 
the symmetry of the problem, Eq. ([5]), and therefore, they could be considered also as the approximate 
solutions of our equations of motion in the limit of large rapidity. The results obtained are impor- 
tant because solutions of RFT based on QCD Pomerons exhibits similar properties and have similar 
structure due the similarities of equations of motion, see [10| . 

The pairs of the functions Eq. ([26]) -Eq. ([27]) and Eq. (f32]) -Eq. ([33]) were known a long time ago, 
see |12j . as "fan" diagrams solutions. They are good approximation for the scattering amplitude in 
the case when the scattering situation is not symmetrical, see [9] [TO]. Therefore, it was interesting 
to check, whether the same functions provide minimal energy for the Hamiltonian also in the case 
of symmetrical scattering with conditions Eq. Q, see "fan" dominance effect discussion in [TO]. Our 
results, thereby, prove that in the semi-classical approximation the "fan" diagrams are the leading 
contribution into the scattering amplitude at the limit of large rapidity also in the case of scattering 
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of symmetrical objects. The stability of these solutions was shown as well by use of the usual stability 
analysis, namely by the Lyapunov functions construction, see Eq. (|54p . 

Considering the symmetrical solution Eq. (|40p -Eq, (j43]) we see that this solution is suppressed in 
the amplitude, see [10] . in comparison with the non-symmetrical, providing non-zero Hamiltonian. 
Therefore, in calculations of different observables of RFT based on QCD Pomerons, O [TU] , in the first 
approximation the leading contribution comes from the "fan" diagrams. The symmetrical solution 
could be considered, therefore, as sub-correction to the amplitude on the semi-classical level. The 
stability of this solution is checked by the construction of Lyapunov function Eq. (|6ip . which was 
obtained in assumption that condition Eq. (|64p is not satisfied. Is this a sign of instability of this 
solution for some values of parameters is not clear. 

Another result of our investigation is a presence of limit cycles solutions which could be demon- 
strated on the base of the system of equations dual to the considered one in the limit of large rapidity. 
In order to understand possible consequences of such solutions on observables of scattering processes 
in RFT at high energies an additional investigation is required. Also, in the future studies, it will 
be important to understand a mechanism of reduction of the limit cycle to the one from the stable 
points. 

A problem of an origin of three solutions during the evolution of the system with rapidity is a 
subject which was not fully considered in the paper. A simplest analysis shows, that new solution are 
arose when the the function p or q begins to be large enough. It happens approximately at rapidity 
y oc — ^ ln(A A) — > 1, but we have no dynamical description of the process similar to usual bifurcation 
picture. It have be noted also, that more complicated boundary conditions lead to more complicated 
picture of possible solutions of the equations of motion, see [33]. A rise of these solutions it is an 
interesting problem which we plan investigate in our further work as well, as the problems mentioned 
above. 
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